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Abstract 

Beyond the rotating-wave approximation, the dynamics of a quantum oscil- 
lator interacting strongly and off-resonantly with a two-level system exhibit 
beatings, whose period equals the revival time of the two-level system. On a 
longer time scale, the quantum oscillator shows collapses, revivals and frac- 
tional revivals, which are encountered in oscillator observables like the mean 
number of oscillator quanta and in the two-level inversion population. Also 
the scattered oscillator field shows doublets with symmetrically displaced 
peaks. 

Key words: two-level system, cavity quantum electrodynamics, adiabatic 
approximation. 
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1. Introduction 

The two-level system interacting with a quantum oscillator has been in- 
tensively studied both theoretically and experimentally. The model is used 
in a wide range of phenomena, especially in atomic physics where it describes 
a two- level atom coupled to a quantized electromagnetic field [H, 0] . It also 
describes the dynamic properties of mixed valence systems with an electronic 
transfer from one part of molecule to another 0, 0, 0| . Such systems are cou- 
pled vibronically and they are in degenerate or quasi-degenerate electronic 
states. The Hamiltonian describes also the ammonia molecule as a pseudo 
Jahn- Teller system [(J. In solid state physics this model describes many 
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situations such as an electronic point defect in a semiconductor [?J or the 
interaction of a dipolar impurity with the crystal lattice (§]. 

More recently the same physical model has been extended to "artificial 
atoms" in a condensed-matter environment with superconductor 



12] and semiconductor systems |l3j. These tiny solid-state devices like flux 
lines threading a superconducting loop, charges in Cooper pair boxes, and 
single-electron spins exhibit quantum-mechanical properties which can be 
manipulated by currents and voltages 14], UM, ll6j. Thus the vacuum Rabi 



splitting has been demonstrated in a Cooper pair box resonantly coupled to 



a cavity mode (lOj, whereas in the dispersive regime, the photon quantum 



number has been probed 12 



The solid-state devices offer wider regimes for the coupling strength be- 
tween the two-level system and a quantum oscillator. Typically, the coupling 
strength in atomic systems is g/u = 10~ 7 — 10~ 6 , where u is the frequency 
of the oscillator, and g is the coupling between the two-level system and the 
oscillator [2]. Similar dipolar coupling in Cooper pair boxes and Josephson 
charge qubits is 3 — 4 order of magnitude larger than the coupling in atomic 
systems |l0|, ll7|. In contrast to the dipolar coupling, the capacitive and 



inductive couplings show even larger coupling strengths 0,0, 19]. 

When the two-level system is conceived as a quantum bit (qubit) which in- 
teracts with a quantum oscillator inside a cavity, the dispersive (off- resonant) 
regime greatly reduces the coupling of the qubit with the environmental noise 



20. 21 



The dispersive regime can be defined as the limit where qubit cavity 
detuning is larger than the coupling (g -C 6). Moreover, recently, it has been 
argued that the regime g/u ~ 1 could be within the experimental realm 



221 ] . Accordingly, the dynamical behaviour of the qubit has been explored 



for the regime v < u (2v is the frequency of the transition between the two 
levels) £ 
coupling (g 



The dynamics of the two-level system show for moderate 
< uj) collapse and revival of the wave-function 22|, whereas 
for very strong coupling (u < g) [23[ the dynamical behaviour is a quan- 



tum version of the Landau-Zener model [24], |25|. Collapses and revivals, the 
adiabatic approximation, and the formal connection between the two-level 
system interacting with a quantum oscillator and the Landau-Zener model 
has been discussed in 



26 



The underlying assumption in deducing the con- 
nection with the Landau-Zener model is a large separation between the two 
displaced harmonic potential curves generated by the interaction between 
the two-level system and the oscillator [23j. Large separation implies large g, 
while for moderate or small g the interference between the dynamics on the 



2 



two displaced harmonic potential curves induces collapses and revivals 22 

Although the dynamics of the two-level subsystem have been studied at 
length in the regime v,g < u, less attention has been paid to the dynamics of 
the quantum oscillator subsystem. On the other hand, system measurements 
can be performed on both subsystems, but probing the oscillator is preferred 
as it is easily accessible via weak measurements 27] . The present study aims 
at the dynamics of the boson field in the regime g < v <C uj and beyond 
the rotating- wave approximation (RWA). The RWA is widely utilized in the 
regime g <C uj, \8\ <C uj, and weak field because the anti-rotating term, which 
is neglected in the RWA, has a quite small contribution to the evolution of 
the system Despite its frequent use in quantum optics, the validity of 
RWA is questionable for applications where strong coupling is present. For 
example, it is known that the energy spectrum of the Rabi Hamiltonian can 
be approximated by its RWA counterpart only for sufficiently small values 
of the coupling strength, and that the width of this range decreases as one 



goes higher up in spectrum [22|, |28|. On the other hand, neglecting terms 
of order g 2 /5 2 within the same RWA, one obtains the shift of the oscillator 
frequency and the quantum counterpart of the ac-Stark Hamiltonian in the 
dispersive regime 2(|. In the present work we will show that the oscillator 
frequency shift, which depends on the qubit state, is obtained without in- 
voking RWA and without satisfying the condition \S\ <C uj. Moreover the 
oscillator frequency shift is asymmetrical , fact that is not found in the RWA 
treatment. 

The regime considered here (g < v <C uj) ensures the adiabatic decoupling 
of the slow motion (the two-level system) from the fast motion (the quantum 
oscillator). The adiabatic approximation shows that the decoupling occurs 
in 2-dimensional subspaces |22f] . We will review the adiabatic approximation 
and we will recast it in a more formal manner that will allows us to use a 
perturbative expansion for the analysis of dynamics. On a medium time scale, 
while the two-level system dynamics shows true collapses and revivals, the 
oscillator dynamics shows beatings. However, on a much longer time scale the 
oscillator itself shows true collapses and revivals accompanied by fractional 
revivals. Fractional revival patterns are present in various observables of both 
the quantum oscillator and the two-level system. A last interesting part is 
the analysis of the scattered oscillator field that shows a doublet which is 
characteristic to dispersive regime. 

The paper is organized as follows. Section I is the introduction. Section II 
presents the Hamiltonian and discusses the adiabatic approximation. Section 
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Ill is dedicated to the analysis of the dynamics. In the last Section the 
conclusions are outlined. 



2. The Hamiltonian and its adiabatic approximation 

2.1. The Hamiltonian 

We consider the interaction of an electronic two-level system with a quan- 
tum oscillator (boson field). In the basis of the two electronic states, |1) and 
1 2), the Hamiltonian can be written as (h = 1): 

p 2 + u 2 q 2 

H = gqa z - va x , (1) 

where o~ z and o~ x are the corresponding Pauli matrices, p and q are the oscilla- 
tor coordinates. The essential parameters are u, g, and v associated with the 
frequency of the oscillator, the electron-oscillator coupling strength, and the 
separation of the electron levels, respectively. This Hamiltonian is isomorphic 
with the simple pseudo Jahn- Teller system 0] or the Rabi Hamiltonian 
[a rotation in a space with e~ l ~ ay or H written in the basis ^ (|1) ± |2))]: 

p 2 + u 2 q 2 

H= - gqa x -va z . (2) 

In this form the parameter v directly shows its role as the splitting of the 
two levels. A simplification of (T5]) is the famous Jaynes-Cummings model 
that has played a very important role in understanding the interaction be- 
tween radiation and matter in the resonant regime [30( ■ With three essential 
parameters, in the dispersive regime, the Hamiltonian ([!]) exhibits a wide 



variety of phenomena like the squeezing of the boson field |3l| , the collapse 



and revival of the wave packet [22J, and the full quantum Landau-Zener ef- 



fect [23]. The simplicity of the Hamiltonian ([T]) led to the conjecture that it 



has an exact solution in terms of known functions [321 . |33| but an analytical 
solution is lacking. Therefore approximations are needed to understand the 
dynamical behaviour of ([1]). Throughout the rest of this paper we are going 
to work with a dimensionless time u t — > t by setting u = 1 without the loss 
of generality. In the next subsections we will present a stationary adiabatic 
approximation of the two-level system strongly coupled with a quantum os- 
cillator. It will be based on the displaced oscillator basis provided by the 
strong coupling and on the relationship between the other two parameters, 
oscillator frequency that is taken to be 1 and the frequency associated with 
the energy separation of the two levels, namely «<1. 
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2.2. Adiabatic approximation to the relative strong coupling (0 <C g < v <C 1) 
The Hamiltonian (HQ) is a two-site realization of the small polaron model 



34j. In the limit v — it can be solved exactly by the Lang-Firsov transfor- 



mation 35]: U = e~ t9peTz . This unitary transformation expresses the Hamil- 
tonian into the basis generated by the displaced oscillators that diagonalizes 
the sector v = of the Hamiltonian pi) . The new Hamiltonian is 

if 1 ( 2 , 2x ^ 2 ( e 2 ^ \ 



Perturbation calculations can be performed as long as v <C 1 [36(. The net 
effect of e ±2%9P on the wave-function is to displace it by the amount ±2g. 
Thus, the effective splitting given by the third term of the Hamiltonian ([3]) 
will be quenched, such that perturbation calculations on the Hamiltonian 
(J2J) can be extended to larger values of v j23|. Going further, we cast the 
Hamiltonian fl3]) in the basis 4g (|1) ± |2)), 

~ _ (p 2 + q 2 - g 2 ) _v_( e 2i9P + e~ 2i9P e~ 2i9P - e 2igp \ 
~ 2 2 \ e 2i9P — e~ 2i9P —e 2i9P — e~ 2i9P ) ' 



So far no approximation has been made. The perturbation calculations that 
consider only the diagonal part of the last term in Eq. (j3j) and the Fock states 
of the oscillator lead to the adiabatic approximation due to the fact that 
the perturbation which is proportional to v is much smaller than the energy 
quanta of the oscillator (v 1). We denote the states of the adiabatic 
approximation as {|1; n) , |2; n)}, where 1 and 2 are the electronic states 
and n is the oscillator quantum number. It is easy to check that the off- 
diagonal part of the last term in (j4j) brings no contribution in the first order 
approximation. Furthermore, expanding the exponentials up to the second 
order in g, we obtain a shift in oscillator frequency. The shift depends on the 
state of the two-level system and it is obtained without using the RWA as it 



has been reported in Ref. [20[. Moreover, in the next section we will show 
that, in contrast to [20], the oscillator frequency shift is uneven. In addition, 
the Hamiltonian (J3J) allows an easy perturbative approach that we will use it 
in the next section. Also, Eq. (j4j) provides a more transparent way to improve 
the adiabatic approximation and to construct a generalized rotating-wave 



approximation (GRWA) as the one that has been recently published 37 



In the formulation given by the Hamiltonian the GRWA is restricted 
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to coupling the states {|2;n) to |l;n + l)}, thus the eigenvalues and the 
eigenvectors of GRWA can be calculated straightforwardly. 

The same adiabatic result is obtained if the matrix elements of Hamilto- 
nian ([1]) are calculated directly in the basis \N±) (Fock states) of the oscilla- 
tors displaced by ±g, where N± is the number of quanta of the state \N±). 
The adiabatic approximation as described above amounts to coupling the 



states \N+) with the states |AL) provided that iV_ = iV+ |22j|. The Hilbert 
space of the adiabatic Hamiltonian is split in a direct sum of 2-dimensional 
sub-spaces {|AL) , |iV+)}. Each sub-space {|iV_) , |iV + )} has the following 
eigenstates and eigenvalues 

\±,N) = -±={\N + )±\N-)), (5) 

E ±tN = ±v (N_ \N + ) + N + (6) 
The energy eigenvalues given by Eq. (jBj) provide a very good approximation 



to the full Hamiltonian Q 221. In the next section we will use the adiabatic 



approximation to explore the dynamics of the quantum oscillator and of the 
two-level system on a time scale of up to O(^r)- 

3. Results and discussions 

The integration of the full dynamics was made with a split-operator 
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method [3JJ. For numerical applications we consider as initial conditions 
a wave packet of the form |1) \a), where |a) is the coherent state [H, IjT" 
which gives the average quantum number of the bosonic part (n) = \a\ 
Collapses and revivals of the two-level wave-function 39j as well as of the 
oscillator wave-function are shown in Fig. [TJ In fact, collapses and revivals 
of the oscillator wave-function are beatings. 

The collapse and revival of the wave-function in the two-level system can 
be understood in terms of the composite Rabi frequency v^ou, i.e. the total 
Rabi frequency is made from contributions of many different Rabi frequencies 
of individual eigenstates of the quantum oscillator (Fig. |2J). When enough 
time has passed, different oscillatory terms get out-of-phase and cancel out 
each other. This is the quantum collapse and it is related to the spread of the 
Rabi frequency distribution that is provided by the wave packet. The revival 
or the restoration of the two-level wave function to nearly its initial value is 
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Figure 1 : Collapse and revivals of two-level (full line) and quantum oscillator (dotted line) 
systems. Only the envelopes are shown. There are two cases: (a) |a| 2 = 30, g — 0.0333, 
v = 0.1; (b)|a| 2 = 100, g = 0.05, v = 0.1. 




Figure 2: Fourier spectrum of (er z ). Low frequency components are characteristic to the 
Fock states of the oscillator, while the high frequency feature shows the coupling to the 
quantum oscillator. 



due to the discrete nature of the Rabi frequency distribution and the revival 
time is related to the difference of the neighbor Rabi frequencies. Revivals 
are purely quantum phenomena, arising from the fact that the quantum 
number distribution of the oscillator is not continuous, whereas the collapse 
is completely classical. The collapse and revival of the two-level system is also 
shown by the adiabatic solution given by Eqs. (jSJ) and (jBJ), whose dynamics 
preserves the general features of the collapse and revival phenomenon as 



just sketched above [22J. The collapse and revival times can be evaluated 
by expanding the off-diagonal exponentials in Eq. ([3]). Thus retaining the 
terms up to second order in p and sandwiching the Hamiltonian ([3]) between 
the adiabatic states {|1; n) , |2; n)} we can estimate the variation of the Rabi 
frequency from level n to level n + 1 

VRabi{n + 1) - v Rabi {n) = Avg 2 . (7) 



Equation (J7J) together with the spread of the wave packet a = v N provide 
the revival time 



and the collapse time 



7T 

tlrev = 2V ( 8 ) 



i cnii - - o- (9) 



4a vg 2 

The dynamics of quantum oscillator, however, show beatings. The beat- 
ings of the quantum oscillator can be understood by analyzing the Fourier 
transform of oscillator coordinate. The spectral features of the oscillator 
coordinate (q) are plotted in Fig. [3j A split-like feature is exhibited in fre- 
quency by the exact dynamics of the quantum oscillator, i.e. the two-level 
system "pulls" the oscillator frequency that now depends on the state of the 
two-level system. 

Quantitative analysis of the splitting can be done by considering the 
diagonal part of the Hamiltonian (j3J) and discarding the constant term —g 2 /2 

~ (p 2 + q 2 ) v ( e 2 ^ + e~ 2 ^ \ 

adiab ~ 2 2 \ - e ^9P _ e -2i<yp J • l iU J 

The expansion of the exponentials in Eq. ffTU]) up to the fourth order in p 
provides the approximate Hamiltonian 
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Figure 3: Spectra of the oscillator coordinate (q) around the unperturbed oscillator fre- 
quency for various coupling strengths g and v = 0.1. They show splittings due to the 
interaction with the two-level system. The finer resolution that is presented in the inset 
indicates the splitting is dependent on the Fock states of the oscillator. 
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H' 



with fin 



j_ 

U 2 2 



\{p 2 + n 2 + q 2 )-v-lvgY 







4^4 



■±(p 2 + tt 2 _q 2 ) + v + pgy 



+0 (p 6 ) 

(11) 



, ~ l=p2vg . The validity of the expansion ffTTj) is bounded 

y/ l±Avg 2 

to maximal quantum number n given by vn 2 g^ < 1, which sets the range of 
validity beyond any practical application as long as g < v <C 1. For instance, 
a typical example v — 0.1, g — 0.1 (largest value of g when i> is set) implies 
n < 316, which is more than enough for a practical application. 

Thus the splitting calculated to the second order in g is Q- — Q + = 4vg 2 , 
i.e. it is just the variation of the Rabi frequency from level n + 1 to level n 
that is given by Eq. (17!) . Therefore the minimum amplitude of the oscillator 
coordinate due to the beatings occurs at the half of the revival time of the 
two-level wave-function, hence the oscillator amplitude may indicate that 
the two-level system is in a pure state [40j. In addition, the fractions in 

Eq. fill I) show the fact that the splitting is uneven. This fact is also found 
in the full numerical calculations shown in Fig. [3J We note here that the 
asymmetrical splitting is not found with a RWA treatment 20 . 



Spreading in quantum number n induces broadening in frequency pulling, 
thus the splitting is well separated only for quite large g (~ 0.1). In the inset 
of Fig. [3] it is presented a finer resolution plot of splitting which shows its 
dependence on the Fock states of the oscillator. The finer resolution will lead 
us to longer time scale dynamics and to true collapses and revivals of the 
oscillator wave-function. The same beatings of the oscillator coordinate can 
be also seen in another regime, v > 1 for a spin coupled to a nanomechanical 
resonator 4l). In that regime, v > 1, the adiabatic motion of the oscillator 
undergoes a frequency change from 1 to yi — g 2 /v 31[ , therefore a classical 
picture of beatings is established with the driving force at the frequency 
of the unperturbed oscillator and the oscillating system at the frequency 
y/l-9 2 /v. 

The analysis can be carried out further by considering the approximate 
wave-function of the Hamiltonian ffTTT) . 



I^i) 



( e -i<p!t J2 C n e - in ( 1 + 2v 9 2 -vg 4 ) t e ~in 2 (-vg 4 ) t |^ \ 
n 

e -if2t ^ C n e - in ( 1 ~ 2v 9 2 + v 9 4 ) t e ~ in2 ( v 9 4 ) t \ n \ 



\ 



J 



(12) 



10 



where C n are the coefficients of the initial wave packet, 0i = — v + |(1 + 
2vg 2 — t>g 4 ), and <\>2 = v + |(1 — 2t>g 2 + t>g 4 ). The validity of the wave- 
function ffl2l) is controlled by the stronger condition vn 2 gH < 1 because the 
energy corrections are included in the phases of the fast oscillating terms. The 



approximations used in [28| and [37[ would describe better the dynamics on 
such large time scales and on wider range of parameters. When using those 
better approximations, however, there are no simple relationships between 
the time scales of the dynamics and the parameters of the system. Thus one 
has to invoke approximations like (fill) and (fl2|) in order to have meaningful 
relationships between the time scales of the dynamics and the parameters of 
the system. It is worth mentioning that a weak coupling limit is used in Ref. 



22j | for the evaluation of the revival time. The wave-function (1121) resembles 
the wave-function of an anharmonic oscillator and it suggests the analysis 
of longer time scale dynamics [42]. The approach of studying collapses and 



revivals as outlined in 42] is, however, suitable for single potential surfaces. 



Once the potentials are coupled, like in the present model, the analysis is 



more elaborate 43]. Thus we pursue an explicit evaluation of the expectation 



values of several observables like 



(<r z ) « 2Re ((V> 2 I Vi» = 2#e ( e -*to-«0* ^ \C n 

\ n 



2 -in(4vg 2 -2vg 4 y 2in 2 vgH 



(13) 

while the expectation value of q up to a residual term g{a z ) has the following 
form 

(<?> - 9 (°z) w (if)i\q \ih) + (H q \ih) = 

73 E v^Q-i^ne-< 1+2V )V^ + ^^Cl +l C n e< 1+2v92 - 2v9i ) t e- 2inv ^ + 

n 
n 

(14) 

Moreover, the difference between the expectation value of the oscillator quan- 
tum number (n) and its time average is 

(n) - (n) t w 2Re (g (ip 2 \ q |^)) = g^Reie'^-^^x 

n * 

(15) 
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Figure 4: The collapse and revival of the oscillator wave-function. The top panel depicts 
the time evolution of {q 2 } that has the revival time halved with respect to (q). The bottom 
panel shows the dynamics of (q) (dashed line), (a z ) (solid line) and (n) — (n) t (dotted line). 
The parameters are: \a\ 2 — 16, g — 0.09, v = 0.1 
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Here (n) t is the time average of (n). For a time scale of order 0(-^^) the 
terms in Eq. (TlBl) rephase, thus one recovers the revival time of (a z ) that 
is given by Eq. (JH]). The examination of Eqs. (jT3"I) . f|T^|) . and (jT5l) reveals a 
longer time scale of order O(-^j^) when their terms rephase again leading to 
the revival of the oscillator wave- function. In Fig. H] we have plotted the time 
evolution of (q 2 ) in the top panel and of (q), (u z ), and (n) — (n) t in the lower 
panel. The new revival time 

7T 

t2rev = 7) (16) 

vg* 

is the revival time of q. The new time t2 re v that is determined by the nearest 
neighbor level spacing is twice the revival time of (q 2 ), which is set by the next 
to nearest neighbor level spacing. Collapses and revivals on a much longer 
time scale have been also predicted in a driven Jaynes-Cummings model 
but the observable subjected to collapses and revivals is (n) instead of 
(q). The case analyzed in the present work shows that, in terms of collapses 
and revivals, (n) behaves similarly to (<J Z ). Moreover t 2 re« is similar to the 



fractional revival time Tp^ identified in the Jaynes-Cummings model [45 
We notice that the time scales defined by the adiabatic approximation are 
different from the time scale defined by RWA and Jaynes-Cummings model. 
In the Jaynes-Cummings model t\ rev and t 2rev depend strongly on the initial 



wave-packet by the mean number of oscillator quanta [45[ . 

A close look in the interval [0, ^rej exposes rich structures of fractional 
revivals in the vecinity of (p/ q)t2 re v (p and q are mutually prime integers) j46[ . 
The analysis of the long time behaviour indicates us that the identification 
of fractional revivals in an oscillator strongly coupled to a two-level system 
can be made by simply analyzing observables like (a z ) or (n) instead of 



calculating the information entropy 47J . Thus in Fig. [5] one can see that at 



the beginning and the end of the interval [0, t2 re v] the collapses and revivals of 
(<j z ) and (n) — (n) t are similar or "in phase" , i.e. they occur at the same time. 
However, in the vicinity of t^revl^ the same collapses and revivals are "out 
of phase": the amplitude of (a z ) is collapsed and the amplitude of (n) — (n) t 
is revived, or conversely the amplitude of (n) — (n) t is collapsed and the 
amplitude of (a z ) is revived. The "in phase" collapses and revivals occur, for 
example, in the vecinity of £2rei>/3, but the revival time is one third of the 
revival time ([8]). The above behaviour is explined by invoking Eqs. fflBl and 
f fTS]) . Thus at t pa t 2r ev/'2 the phases brought in by the terms e 2mvg t e 2m vg * 
in Eq. ( 1131) are 0. In turn the phases brought in by the terms e 2m vg 1 in 
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Figure 5: (Color online) Fractional revivals that can be observed in (a z ) and (n) — (n) t . 
(a) and (b) at the beginning and the end of the interval [0,t2 ret ,], respectively; (c) around 
^2ret>/2; (d) around t2reu/3; (e) and (f) in the vicinity of ^2re«/4 and 3<2ren/4, respectively. 
The thin black lines denote the dynamics of (a z ) and the thick red lines depict the dynamics 
of (n) — (n)t- The parameters are the same as the ones used in Fig. S) 
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Figure 6: (Color online) Full dynamics of (a z ) and (q) in the vicinity of t^rev/^- The full 
line depicts the dynamics of (q) and the dotted line denotes the dynamics of (c z ). 

Eq. f[T5]) are ir. Consequently the collapses and revivals of (a z ) are "out of 
phase" with respect to collapses and revivals of (n) — {n)t- Similar behaviour 
is shown at t pa t2rei>/4 and t ~ 3t 2rp „ /4, while at t ~ ^2re^/3 the revival time 
of (a z ) and (n) - (n) t is h rev /3 [46J. 

Another interesting aspect is the frequency of (q) within collapse region 
as it can be seen in Fig. EJ the frequency of (q) equals the frequency of 
the two-level system. Meanwhile in the revival regions the frequency of (q) 
is close to the frequency of the unperturbed oscillator. The inspection of 
Eq. (fbil) tells us that, in the collapse region, (q) is mainly the term g(a z ), 
hence the dynamics of (q) and (<j z ) are similar in the collapse region of the 
oscillator. 

The approximate wave-function f[T2|) allows us to calculate also the scat- 
tered field of the oscillator by the two-level system. We define the scattered 
field as the Fourier transform of the first order correlation function: 

(oo oo \ 

j dt j dt' (at (t) a z {t) a (f) a z (£')) e^^e (t -t')\. (17) 

Direct calculations lead us to 

I scatt oc n \ C nf 6 " 1 =F 2v =F V ± 4 ™V T 2n 2 vg 4 ). (18) 
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Unlike the spectrum of (g), the spectrum of the scattered field shows two 
lobs symmetrically displaced by the low frequency components of Rabi fre- 
quencies with respect to the frequency of the unperturbed oscillator. The 
spectrum of scattered field is the counterpart to the fluorescence spectrum 
obtained by strong field pumping of a two- level system near the resonance [l[ . 
Thus, experimentally the scattered field can be probed in the same way like 
the probing of fluorescence spectrum Fluorescence spectrum has three 
peaks, one dominant central peak at the frequency of the pump field and two 
side peaks symmetrically displaced by Rabi frequency from the central peak 



48( | . The spectrum given by Eq. (fl8l) . however, has two lobs that are charac- 
teristic to a dispersive regime. Moreover, the scattered spectrum shows also 
information about the occupation of oscillator Fock states. 

4. Conclusions 

We have studied the dynamics of a quantum oscillator interacting strongly 
with a two-level system in a dispersive regime without invoking the rotating- 
wave approximation (RWA). First we recast the adiabatic approximation 
by applying two successive unitary transformations. We also pointed out 
the frame toward the generalized rotating-wave approximation (GRWA) that 



improves the adiabatic approximation [37 . 

Collapses and revivals in the two-level system are accompanied by beat- 
ings that are exhibited by the dynamics of the quantum oscillator. The 
beatings are related to collapses and revivals of the wave function in two- 
level system. However, on the much longer time scale the quantum oscillator 
shows true collapses and revivals that occur in the dynamics of the oscilla- 
tor coordinate. Between the new collapses and revivals there are fractional 
revivals whose signature is present in both the two-level system and the 
quantum oscillator. 

Finally, we calculated the scattered oscillator field off the two-level sys- 
tem. Due to the dispersive regime, the scattered field spectrum shows a 
doublet symmetrically displaced about the frequency of the unperturbed os- 
cillator by the amount of the low frequency components of Rabi frequencies. 
The spectrum can be used not only to assess those Rabi frequencies but also 
to probe the distribution of the oscillator quantum number. 
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